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We present a detailed numerical study of the effect of a disordered potential on a confined one- 
dimensional Bose-Einstein condensate, in the framework of a mean-field description. For repulsive 
interactions, we consider the Thomas-Fermi and Gaussian limits and for attractive interactions the 
behavior of soliton solutions. We find that the disorder average spatial extension of the stationary 
density profile decreases with an increasing strength of the disordered potential both for repulsive 
and attractive interactions among bosons. In the Thomas Fermi limit, the suppression of transport is 
accompanied by a strong localization of the bosons around the state fc = in momentum space. The 
time dependent density profiles differ considerably in the cases we have considered. For attractive 
Bose-Einstein condensates, a bright soliton exists with an overall unchanged shape, but a disorder 
dependent width. For weak disorder, the soliton moves on and for a stronger disorder, it bounces 
back and forth between high potential barriers. 

PACS numbers: 03.75.Kk,03.75.-b,42.25.Dd 



I. INTRODUCTION 

The spatial behavior of a wave submitted to a strong 
enough random potential remains one of the major and 
still unsolved issues in physics. It is an ubiquitous prob- 
lem that shows up in almost all fields ranging from astro- 
physics to atomic physics. The interference induced spa- 
tial localization of a wave due to random multiple scat- 
tering has been predicted and named after Andersoni. 
The Anderson localization problem despite its relatively 
easy formulation has not yet been solved analytically and 
still rises a lot of interest. Strong Anderson localization 
of waves has been observed in various systems of low 
spatial dimensionality where the effect of disorder is ex- 
pected to be the strongesliSi^. Above two dimensions, 
a phase transition is expected to take place between a 
delocalized phase that corresponds to spatially extended 
solutions of the wave equation and a localized phase that 
corresponds to spatially localized solutions. The descrip- 
tion of this transition is mainly based on an elegant scal- 
ing formulation proposed by Anderson and coworkers^ 
Due to its indisputable importance, the localization of 
light is a hotly debated but still unsolved problemSi^. 
The weak localization regime, a precursor of Anderson 
localization for weak disorder, has been studied in detail 
both theoretically and experimentally for a large variety 
of waves and types of disorderiSiiSiiiii^. 

In contrast, relatively little attention has been paid 
to the extension of Anderson localization to a non-linear 
medium. Though analyticaliiMii^ as well as numeri- 
cal work have been done to address this issue, no clear- 
cut answers have been obtained to ascertain how local- 
ization is affected by the presence of a non-linear term 
in a Schrodinger type wave equation. This is the pur- 
pose of this paper to address this issue in the context 



of the behavior of a one-dimensional Bose-Einstein con- 
densate (BEG) in the presence of a disordered optical 
potential, since it has raised recently a great deal of 
interesti6,i7d8,i9i20,2ii22,23,24,25,26,27,28^ Transport of a 

magnetically trapped BEG above a corrugated microship 
has been theoretically studied recentlji^S. The possibil- 
ity of tuning random on-site interaction has also been 
considcredi^. Using Feshbach resonances, it is possible 
to switch off the interaction among bosons which will 
then be allowed to propagate through a set of static im- 
purities created by other species of atom. This may lead 
to an experimental realization of the Anderson localiza- 
tion transition. The corresponding theoretical model has 
been proposed and analyzed^^*2^ for one-dimensional sys- 
tems, i.e. in the absence of transition. The other issue is 
to understand the interplay of interaction induced non- 
linearity and disorder on the Bose-Einstein condensate. 
One-dimensional systems are especially interesting since 
the effect of disorder is the strongest and such systems 
are experimentally realizable. Experiments in this direc- 
tion have been performed recentlj*iSiiLiS which show a 
suppression of the expansion of the BEG cloud once it is 
released from the trap. 

In this paper we present a numerical study of the ef- 
fect of a disordered potential on one-dimensional conden- 
sates with either attractive or repulsive interaction in the 
framework of the mean-field approximation and compare 
between these two cases. Studies of the propagation of 
a quasi one-dimensional BEG in a disordered potential 
have been carried out mostly in the repulsive Thomas- 
Fermi li];nitiSiiLi2i2242&. We also consider this limit and 
we find numerical evidence that the suppression of the 
BEG expansion after the release from the trap, is due to 
localization in momentum space around the state fc = 0, 
with becomes stronger for an increasing strength of dis- 
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order. This suggests that the momentum spectroseopy 
of disordered quasi one-dimensional BEC may give im- 
portant information about its transport properties. In 
addition, we consider the Gaussian Umit of a strong con- 
finement and the bright soUton solution for an attractive 
interaction. We also propose a model for the disorder 
where both the strength and the harmonic content can 
be independently varied. 



The comparison between the different cases is moti- 
vated by the fact that depending on the strength and 
sign of the effective interaction among bosons in an ef- 
fectively one-dimensional BEC, various types of scenarios 
may be realized. The interplay between these different 
types of interaction and disorder should lead to different 
types of stationary as well as time-dependent behavior of 
the density profile. We consider three such regimes that 
cover both the repulsive and attractive interaction and 
where the system can indeed be well described within 
the mean-field approximation. The corresponding mean- 
field is given by the Gross-Pitaevskii equation with mod- 
ified coupling constant^ (in comparison to the three di- 
mensional case) and it takes the form of a non-linear 
Schrodinger equation. Its solutions in the absence of dis- 
order have been thoroughly studied2^i2S42L2&. We employ 
a numerical scheme^MSi^L^ which has been recently de- 
veloped to study stationary solutions of this non-linear 
Schrodinger equation in the presence of a disordered po- 
tential. The scheme is based on a rapidly converging 
spectral method. Then we look at the time evolution of 
the stationary profile after switching off the trap poten- 
tial. Subsequently, we analyze our solutions and compare 
them to those obtained in the absence of disorder. Our 
study unveils an interesting picture for the interplay be- 
tween the nature and strength of interaction and a ran- 
dom potential. 

The organization of the paper is as follows. In sec- 
tion II we briefly review the stationary density profiles 
of an effectively one-dimensional BEC in the absence of 
disorder and in the mean field regime. Then, in section 
III, we introduce our numerical scheme and we define our 
model of disorder on such one-dimensional condensates. 
In section IV, we present our numerical results for the 
Thomas- Fermi limit. In section IV- A we compare our 
results with recent works on this subjectiSiiLiSiSfliSSiSS. 
In section V, the effect of disorder in the confinement 
dominated Gaussian regime is discussed. Both sections 
pertain to the situation of repulsive interaction among 
bosons. In section VI, we discuss the effect of disorder 
on a bright solitonic condensate which corresponds to 
an attractive effective interaction. In the last section we 
summarize and present the general conclusions derived 
from our results. 



II. STATIONARY SOLUTIONS IN THE 
ABSENCE OF DISORDER 

A. One-dimensional repulsive Bose-Einstein 
condensate in a trap 

We review briefly the mean field description of a quasi 
one-dimensional Bose gas with short range repulsive in- 
teraction, in a cylindrical harmonic trap along the z-axis, 
and in the absence of disorder. Details are given in 
references^i^. The Gross-Pitaevskii equation provides 
a mean field description of the three dimensional inter- 
acting gas and it is given by 
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where Wz and lo±_ are respectively the harmonic trap fre- 
quencies along the z-axis and along the radial direction; 

az = and a± = \J'^^^ are the corresponding har- 

monic oscillator length scales. The interaction is charac- 
terized by the s-wave scattering length a, that is positive 
for a repulsive interaction. For tight trapping conditions 
[lOz ^ all atoms are in the ground state of the har- 
monic trap in the radial direction and the condensate is 
effectively one-dimensional. Nevertheless, for a± > a, 
the effective coupling constant along the z-direction is 
still characterized by a and it is given by gid — 2afiw_L— . 
The corresponding mean-field behavior is governed by 
the Gross-Pitaevskii equation. 
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where is the condensate wavefunction along the z-axis. 
We look for stationary solutions of the form '^{z,t) = 
(j){z) exp(— i/it) where jx is the chemical potential. The 
corresponding one-dimensional density is pid = 
The interaction strength may be expressed in terms of 
the dimensionless coupling constant 7, 
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which is the ratio of the mean-field interaction energy 
density to the kinetic energy density. For 7 <C 1, the 
gas is weakly interacting and, in contrast to higher space 
dimensionalities, in one-dimension the gas can be made 
strongly interacting by lowering its density. For larger 
values of the interaction strength 7, the Gross-Pitaevskii 
equation (0) does not provide anymore a correct descrip- 
tion, the gas enters into the Tonks- Girardeau regime^^ 
and behaves like free fermions. 

Starting from 0, a dimensionless form can be 
achieved that is given by 
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where use has been made of rescaled length and time, 
z ^ —, t ^ Hf-t and 4" ^/a7^- The dimensionless 
parameter aid, or equivalently the coherence length is 
defined by 
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and it accounts for both interaction and confinement. By 
rescaling the chemical potential, fi — + jjf^, we obtain 
for the time-independent Gross-Pitaevskii equation the 
expression 



/i0 + dl(t> - - 2aid|0p0 = 
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Henceforth we shall express our results in terms of these 
dimensionless quantities unless otherwise specified. We 
mention now two limiting regimes of interest that can be 
described by Eq.lQ. 



1. Thomas-Fermi limit 

For a chemical potential jl larger than the level spacing, 
namely for jl ^ hw^ (i-e. in dimensionless units /i 3> 1), 
the gas is in the Thomas-Fermi regime. Thus the kinetic 
energy term becomes negligible. We denote by ptf and 
jiTF the corresponding condensate density and chemical 
potential. We have 
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The number of bosons is given by A'' = J^^'' dzpxp, 
where Ltf — \/I^tf is the Thomas- Fermi length. Elim- 
inating Ltf, we obtain, 
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2. Gaussian limit 

The other limit jl <^ fu-Oz, corresponds to a regime 
where the single particle energy spacing is larger than 
the interaction energy so that the gas behaves like N 
bosons in a harmonic trap potential. Thus, wc have an 
ideal gas condensate with a Gaussian density profile. 

As we shall see later, the effect of a disordered potential 
on the condensate dynamics for both limiting cases is 
significantly different. 



B. One-dimensional attractive Bose-Einstein 
condensate 

We also consider the case where the s-wave scatter- 
ing length a is negative. The effective interaction among 



bosons is thus attractive. This situation can also be de- 
scribed by means of Eas. l|2l4|l . In the absence of con- 
finement and when aid = —1—, Eq.Q) admits a moving 
bright soliton solution of the form 
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where Vs and p > Q are respectively the velocity and the 
chemical potential of the soliton solution and (zq, </'o) 
refer to the translational and global phase invariance of 
Eq.Q. In particular, if 14 = and choosing for simpHc- 

ity the gauge zq = (j)Q = 0, then 'J'(z,t) = (f)[z) e}cp[ipt) 
with 

(j){z) = sech(^ z) , (10) 
time-independent nonlinear 



which satisfies the 
Schrodinger equation 



-p<P{z) + dl<P{z) + 2\<P\''^ = 0. 
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The chemical potential p is proportional to the square of 
the inverse width of the soliton. Such a soliton has been 
experimentally observed2& and theoretically studied^ for 
cold atomic gases. 



III. NUMERICAL METHOD FOR DISORDER 
AND NON-LINEARITY 

A. Spectral method 

We start by considering the dimensionless time- 
independent Gross-Pitaevskii equation 
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in the presence of a disorder potential Vd{z). Upon dis- 
cretization, this potential is defined at each site of a lat- 
tice and it is given by the product of a constant strength 
Vm times a random number u which is uniformly dis- 
tributed between and 1. Using a Gaussian approxi- 
mation with mean a (the lattice spacing), the disorder 
potential can be written as a continuous function 



with 



V{z — z') ~ lim Vm exp 



Vd{z,u)^ujV{z) , 
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A disorder potential generated in this way varies rapidly 
over a length scale of the order of a lattice spacing. We 
wish however to use a smoother potential more appropri- 
ate for the description of typical disorders generated in 
experimentsiSiiS. To that purpose, we consider the dis- 
crete random variable u) defined at each lattice site and 
we discard from its Fourier spectrum all wavenumbers 



that are above a given cutoff kc = 2tt/Xc. The inverse 
Fourier transform uj{Xc) = provides a random poten- 
tial that varies on length scales larger or equal to Ac and 
which can be formally written as 



By Fourier transforming Ea. (|16|l we obtain 
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where M is a large enough number. The new random 
variable ll!'^(Xc,z) thus generated is different from uj. 
Whereas the average value of uj is, by definition, equals to 
1/2, we obtain, for example, that for kc = 6, the average 
value of iu'^ is about 2 x 10^^. Typical examples of such 
slowly varying potentials obtained by changing Ac are 
given in section V in Fig|5{a,c,e). The disorder potential 
= Vr,iU!''{Xc, z) that we consider is thus characterized 
by two quantities: its strength Vm and the scale Ac of its 
spatial variations. Ea. H12|) rewrites 

lJL(t>u + dl(p^ - z'^cj)^ - Kn'^''(Ac, z)(l)^ - 2aid\<i)uj\'^ 4>u . 

(16) 

The local density for a given realization of disorder is 
Ptjj{z) = |0ij(z)p and the number N of bosons is deter- 
mined by the condition N — J dz Pu{z). By direct in- 
spection of the different terms that show up in Ea. (|16|l . 
we see that disorder effects are obtained either by com- 
paring them to interactions, i.e., by comparing the dis- 
order length scale Ac to the coherence length ^ defined in 
(|5|). If the ratio Ac/^ is small, disorder is strongly varying 
spatially and its effect overcomes that of interactions. We 
also compare the effective disorder strength Vm^'^ to the 
chemical potential This can be achieved by defining 
the local dimensionlcss random variable 



(17) 



We will consider its average over configurations denoted 
by (s). The parameter s allows to compare between the 
chemical potential and heights of barriers of the disor- 
der potential. This parameter, as we shall see, plays also 
an important role in the study of the time evolution of 
the density once the trapping potential is released. It 
controls the spatial extension of the cloud as a function 
of time. Finally, we consider boundary conditions for 
Ea. (|16|l obtained by demanding that for a given realiza- 
tion of disorder, (j)uj{z) vanishes for \z\ — > -foo. 

We now turn to the description of the numerical 
method used to solve Ea. H16|) . The fact that it is ran- 
dom, makes it very challenging for conventional numer- 
ical schemes to be implemented. The numerical scheme 
we use here is based on the spectral renormalization 
method that has been recently suggested by Ablowitz 
and Musslimani^ (see alsc42i^) as a generalization of 
the so-called Petviashvili method^. Spectral renormal- 
ization is particularly suitable for this type of problems 
for its ease to handle randomness. To this end, for a fixed 
realization, we define the Fourier transform 



(j^uiik) = T[<j)u,{z)] = / dz(t)^{z)e 



-ikz 
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2aidJ'[\ct>^\''(l>u,] + J^jz^u^jz)] + T[V%z,u;)(j>^ 
^ — fc2 



fj, — K~ 

(19) 

In general, the solution of this equation is obtained by 
a relaxation method or successive approximation tech- 
nique where one gives an initial guess and iterates until 
convergence is achieved. However, this relaxation pro- 
cess is unlikely to converge. To prevent this problem, 
we introduce a new field variable tpuiiz) using a scaling 
parameter , 

(pujiz) =Pujipuj{z) , 4'u^ik) = Pu-ip^{k) . (20) 

Substituting into Eq. H19(l and by adding and subtracting 
the term r(j)^{k) (with r > 0) to avoid division by zero, 
we obtain the following scheme 
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where pi™'' are given by the following consistency condi- 
tion 



(22) 

where the inner product in Fourier space is defined by 
(/,.9> = / fgdk . 



B. Time dependent evolution 

To describe the time evolution of the stationary so- 
lutions, we use a time splitting Fourier spectral method 
that has been described in detail^. We describe it briefly 
with a comment on its limitation. 

After switching off the trap, the time evolution is gov- 
erned by the equation 

idt-^Uz,t) = -d^,^Uz,t) + V„,uj'{Xc,z)-^Uz,t) 

+2aid|^'^(z,i)|2*^(z,f) , (23) 

with ^'j^(z,0) = (j)uj{z). Equation H23|l is solved in two 
distinct steps. We solve first 

idt^Uz,t) = ^d^,^Uz,t), (24) 

for a time step of length At and then, 

idt'fUz,t) = V^uj^Xc,z)-f^{z,t)+2aid\'f^{z,t)\^'if^{z,t), 

(25) 
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for the same time step. The first of these two equa- 
tions, (|24|l . is discretized in space by the Fourier spectral 
method and time integrated. Tire solution is then used 
as the initial condition for the second equation H25|) . The 
commutator between the two parts of the Hamiltonian 
that appears in the right hand side of H24|l and (|25|l is 
disregarded in this process. The resulting error is signifi- 
cant if this commutator is large compared to other terms 
in the equation. This is the case if the disordered poten- 
tial strongly fluctuates (which is not considered in the 
present numerical work). Notice that, by definition, this 
method ensures the conservation of the total number of 
particles. 



IV. THOMAS-FERMI LIMIT 
A. Stationary solutions 

Stationary solutions to the Gross-Pitaevskii equation 
(ini in the presence of a disordered potential and in the 
Thomas- Fermi limit can be obtained by iterating Eqs. 
H21(l and H22() . Then, we compare these solutions with 
those obtained by directly considering the Thomas- Fermi 
approximation in the presence of disorder. This compar- 
ison is displayed on Fig^ 




FIG. 1: (a) Behavior of the condensate density \(j>uj\'^ obtained 
from the Thomas-Fermi approximation and from the Gross- 
Pitaevskii equation for a spatially slowly varying disordered 
potential characterized by Ac ~ 12^. This corresponds to a 
weak disorder as compared to interactions, (b) Same plot but 
for a spatially rapidly varying disorder such that Ac ~ 2^. The 
average disorder is kept much below the chemical potential 
(yL — 30) in both cases, so that (s) <C 1. We have taken 
Old — 1 and a number N of bosons equal to 80. 

Generalizing the Thomas-Fcrmi approximation ^ so 



as to include the disorder V^, we obtain for the corre- 
sponding density the expression. 
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It is thus expected that this density presents local max- 
ima and minima that follow the corresponding ones of the 
disordered potential. This trend is indeed very apparent 
in Fig^ In Yig.\.(a), we consider a very smooth disorder 
such that Ac is much larger than the coherence length ^ 
and we observe that apart from little deviations, the den- 
sities obtained from the Gross-Pitaevskii equation and 
from the Thomas-Fermi approximation match almost ex- 
actly as expected. Such a smooth disorder corresponds 
to the typical situations encountered in the experiments 
performed at OrsajiiS, at LENSii and Hannover—. In 
¥ig.l.(h), the disorder is stronger i.e. that it fluctu- 
ates on a smaller length scale Ac comparable to ^, thus 
leading to more local minima and maxima of the dis- 
ordered potential within the size of the cloud. The 
Gross-Pitaevskii density, obtained by solving (|16|l . de- 
viates from the Thomas-Fermi density at these extremal 
points. Moreover, we observe that the magnitude of those 
deviations gets larger when Ac gets smaller, i.e., for larger 
spatial variations of disorder. This behavior can be un- 
derstood by considering the following expression for the 
density p^^ 
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which follows straightforwardly from Ea. ljlGII . In this ex- 
pression, the second term in the r.h.s., also known as 
quantum pressure term, is a correction to the Thomas- 
Fermi density whose origin is the zero point motion of the 
bosons in the condensate. This correction is proportional 
to the ratio (^/Ac)^. It becomes larger for a decreasing Ac, 
namely for a relatively larger effect of interactions driven 
by ^. Thus a stronger disorder introduces more apprecia- 
ble zero point motion of the bosons so as to reduce the 
interaction energy cost. In other words, the behavior of 
the static Thomas-Fermi condensate in a random poten- 
tial is such that the disorder potential gets screened by 
the repulsive interactionSS*^. 

Another feature of disorder is the spatial extension of 
the cloud defined, for a given disorder configuration, by 



(28) 



where we have characterized the spatial distribution of 
the cloud by its moments. 
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In Fig[21 wc have plotted the configuration average {LJ) 
of the spatial extension as a function of the average 
strength ^(s) (see Ea tTTj) . The average spatial extension 
of the cloud, in the Thomas-Fermi limit is a decreasing 
function of the ratio Ac/^, i.e., it decreases when inter- 
actions are getting larger than the spatial variation of 
disorder. We shall see that this behavior holds true even 
beyond the Thomas-Fermi approximation. 



comes spatially localized. Then, the average parameter 
(s) is not anymore relevant since it may be smaller than 
unity although some barriers may be larger than /i. We 
thus need to characterize the disorder by means of higher 
moments. For a smooth enough probability distribution 
of the random variable KnW'^, which is the case we con- 
sider, it is enough to consider the variance Slo'^ defined 



by Slo" = (((w^) ) - K)- 
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and the parameter 
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FIG. 2: Plot of the width (Lu,) averaged over disorder, as a 
function of the average strength fJ,{s). (L^) is expressed in 
units of the corresponding extension in the absence of disor- 
der. For a given value of Xc, we average over 200 realizations 
of the disordered potential. Such calculations have been done 
for five different values of Ac, the last one being less than 
the coherence length ^ for which the validity of the mean-field 
theory is questionable. We have taken /i = 30 and aid = 1 



B. Time evolution 

We study now the time evolution of the previous sta- 
tionary solutions while switching off the trapping po- 
tential, but keeping the disordered potential. We then 
compare them with recent experiments and numerical 
calculationsiSiiii2i2£. In the experimentsi^iSa, the BEC 
was prepared within the trapping and random potentials, 
but its expansion has been studied while switching off 
these potentials. This led to the observation of sharp 
fringes in the resulting density due to interference be- 
tween different parts of the condensate. These conditions 
therefore differ from the case we consider. 

We recall that our disorder is characterized by its 
strength s in units of the chemical potential fi and by 
the length scale Ac of its spatial variations. The latter 
quantity is analogous to the correlation length of disorder 
defined iniS. It is important to stress that in the Thomas- 
Fermi regime, the time evolution is very sensitive to the 
existence of potential barriers of height larger than the 
chemical potential /i. If such a barrier exists, say at a 
point zq, then we observe that the density Pujiz) vanishes 
for z > Zq at any subsequent times so that the cloud be- 
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which sets the width of the distribution of potential bar- 
riers. In some of the cases we consider, the peak height 
of the disordered potential is twice as high as Kn^^j^. A 
specific feature of one-dimensional disorder is that it is 
always very strong in contrast to higher dimensional sys- 
tems for which the cloud may always find a way to avoid 
large potential barriers thus making effects of disorder 
comparatively weaker. 

In Fig|3|we have plotted the density profile Pu{T) after 
a time T for different spatial variations of disorder. For 
instance, Yig&.'i(a) and (b) compare cases with different 
values of Ac but keeping the disorder strength 5s <C 1, 
almost unchanged. In contrast, Figl^J^aJ and (c) display 
time evolutions for two disorder potentials having the 
same value of Ac, but different strengths 5s. In FigUfc), 
one of the potential well has a height almost equal to 
/i. A first general observation is that for smaller values 
of Ac, namely for stronger spatial fiuctuations, the spa- 
tial expansion of the cloud is more inhibited, so that the 
main part of the cloud remains localized in finite regions 
that depend on the local landscape of the disordered po- 
tential. This trend is clearly apparent in Fig|2| On the 
other hand, for small values of 5s, small amplitude den- 
sity fiuctuations extend far apart from the initial point 
whereas for large values of 5s, density fluctuations do not 
extend beyond the closest barrier of height larger than /i. 

It has been pointed out that, during the time evo- 
lution of the Thomas-Fermi cloud, its center and its edge 
behave in a different way. After the trap potential is re- 
leased, the density peak at the center, that corresponds to 
the highest value of the stationary density, gets lowered 
at an initial stage of the expansion. The interaction en- 
ergy being larger than the kinetic energy, the density pro- 
file near the center still closely follows a Thomas-Fermi 
shape, but with a reduced chemical potential. The spa- 
tial variation of density fluctuations corresponds approx- 
imately to that of the disordered potential (that is of 
the order of Ac). At the edges of the cloud, the density 
is lower so that the kinetic energy term takes over the 
interaction term and it is almost equal to the chemical 
potential p, of the condensate a,t t — Q. Thus, the charac- 
teristic scale of spatial variations of density fluctuations 
at the edges of the cloud is the coherence length ^ which 
is smaller than Ac. This is displayed in Fig0] which de- 
picts the time evolution at the center {(a) and (b)) and 
at the edge {(c) and (d)). In Fig^ (a), the center of the 
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FIG. 3: (a) Plot of the density profile — \cI>lu{z)\'^ for two 
different times, T = and T — 25/ujz and for a slowly vary- 
ing disorder such that Ac ~ 12^. Both the average strength of 
the potential (s) = 0.104, and 5s — 0.078 are kept much below 
unity and ^ = 30. The variables plotted in the vertical axis are 
indicated on the right hand side of the upper figure. The den- 
sity profile in the absence of disorder at T = 25/uJz is shown 
in all three figures. The cloud extends almost uniformly. The 
horizontal axis label is indicated in the first figure only, (b) 
Time evolution of the density for a stronger spatial variation 
of disorder such that Ac ~ 2^. (s) = 0.095 and 5s — 0.076 
are kept well below unity (jj, = 30^. (c)Time evolution of the 
density for a stronger disorder characterized by (s) — 0.23, 
5s = 0.185 and Ac = 12^. aid = 1 for all three plots 



time evolved cloud follows the potential landscape and 
varies on a much larger length scale than the edge of the 
cloud. The other limit shown in Fig0|(('&j and (d)), dis- 
plays relatively less difference between spatial variations 
of density fluctuations at the center and at the edges of 
the cloud, since Ac — ^. 

In addition. Figs 01 describe how the matter wave be- 
haves close to a single potential barrier, at the center and 
at the edge of the cloud. The shape of a typical potential 
is controlled by changing Ac- In Figs^ (a) and (b), it is 
shown how the central cloud becomes localized due to the 
presence of a potential barrier. The density modulation is 
driven by the the local potential landscape, rather than 
any interference effect. It has been pointed out iniLSS 
that the height of a single defect should vary like the en- 
ergy E of the incoming wavepaeket over a distance short 
compared to its de Broglie wavelength in order to allow 
for quantum effects to dominate and eventually lead to 
Anderson localization. The potential used in our com- 
putation, generally does not satisfy this criterion. To 
satisfy it, one needs a disorder with higher Ss and lower 
Ac. However under such conditions, the mean-field Gross- 
Pitaevskii approximation is questionable and the use of 
a discrete non-linear Schrodinger equation will be more 
appropriate. 



(a) center (b) p (t> 




z 4 4 




58 2 62 38 42 



FIG. 4: (a) and (c) Time evolution of the density at the center 
and at the edge for a strong disorder ( (s) — 0.23, 5s — 0.185 
and Ac ~ 12^j identical to the one used in FialWc). (b) and 
(d): Time evolution of the density at the center and at the 
edge for a disorder characterized by (s) = 0.095, 5s — 0.076 
and Ac = 2^, identical to the one used in Fig\^b). ^ = 30 and 
aid = 1 for all figures. The horizontal, vertical and color-axis 
are the same for all three plots and are shown in alternative 
pair of figures. The black line in each figure shows the cor- 
responding disordered potential. The potential is rescaled and 
its origin is shifted by the same amount in all figures to fit it 
in the size. 



We have studied in Fig|Sl the time evolution of the 
cloud density in momentum space and compare it to the 
cases without disorder and in the presence of an optical 
lattice. To make the comparison easier, we have used the 
amplitude of the optical lattice potential corresponding 
to that, < s > -|-(5s, considered in Fig[Sl (a). Fig|31 (a) 
shows a strong localization in fc-space for high 6s. This 
may be compared to the situation of FiglSJ^dj (optical 
lattice). This strong localization occurs aroimd the k ~ 
state. On the other hand when disorder fluctuates on a 
shorter scale Ac (with a smaller 5s), a, significant fraction 
of the density still occupies higher momentum states and 
the corresponding localization in momentum space is less 
pronounced. Thus, an experimental measurement of the 
momentum spectrum^ of a quasi one-dimensional BEG 
in a disordered waveguide can shed light on the nature 
of localization of the cloud. 

After studying the time evolution of the density, we 
now discuss the time evolution of other properties of the 
cloud that characterize the suppression of its expansion. 
In FigJSJa), the spatial extension for a given configu- 
ration, is plotted as a function of the dimensionless time 
uj^t. We observe that L^{t) saturates with time to a 
value which depends on the average strength (s) of the 
disorder. 

In order to characterize this saturation, we define the 
ratio, denoted by TZ, between the average kinetic and 
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(a) 1^ (k,t)| 




FIG. 5: Density evolution in the k — space. Horizontal and 
vertical axis labels are identical for all the plots and are shown 
in (a). The color-axis label is indicated in (a) and (c). The 
four cases correspond to (a) Strong disorder as defined in 
Fia\^c). (b) Weaker disorder fluctuating on a smaller length 
scale as defined in Fig^(b). (c) No disorder (d) Optical lat- 
tice: ((s) + 5s) sin . The values of{s), 5s and Xc are those 
used in Fig\^c). 



interaction energies of the cloud, defined by 



7^ = 2C 



Sdz\<t>X 





(31) 



no disorder 

[<s>, A, , 5s] 
0.104,1*^^,0.078 

0.23, 12^,0.185 

0.095, 2^,0.076 



FIG. 6: (a) Time evolution of the spatial extension Lu]{t) of 
the cloud as defined in i2St for a given configuration of the 
disordered potential but for different strengths V,n and length 
scales Ac. (s), Ac and Ss are indicated against each plot. L^{t) 
is expressed in units of its value at t = 0. (b) Corresponding 
time evolution of the ratio TZ defined in 



V. GAUSSIAN LIMIT 

In this section we study effects of disorder on bosons 
that are condensed in the ground state of a harmonic 
oscillator potential. In that case, solutions of the Gross- 
Pitaevskii equation without disordered potential, are dif- 
ferent from those observed in the Thomas-Fermi limit, 
and are given by Gaussian profiles centered at the origin. 
When the trap potential is released, corresponding time- 
dependent solutions remain Gaussian but with a larger 
width, a standard result from quantum mechanics (see 
FigHa)). 



In the stationary Thomas-Fermi approximation, the ki- 
netic energy is almost negligible as compared to the in- 
teraction term, i.e., TZ ~ 0. As the cloud expands, the 
interaction energy gets gradually converted into kinetic 
energy and this ratio increases. Finally, it saturates when 
almost all the interaction energy is converted into kinetic 
energy. This shows up in Fig|^6). For a larger disorder, 
this increase of the ratio saturates more rapidly and the 
slope of TZ, which indicates how fast the interaction en- 
ergy is converted into kinetic energy, decreases. Partic- 
ularly the lowest plot corresponding to a large disorder, 
shows a rapid saturation of TZ due to strong localization 
in momentum space. Since the edge of the cloud involves 
mostly kinetic energy, the behavior of TZ is dominated by 
the expansion of the central region. When the expansion 
is stopped by a potential barrier, the corresponding loss 
in kinetic energy is proportional to the height of the po- 
tential barrier. This explains the oscillations of TZ that 
appear in the presence of disorder. 



A. Stationary solutions 

Like for the Thomas-Fermi regime, stationary solutions 
of the Gross-Pitaevski equation H16|l in the presence of 
both trapping and disorder are characterized by the av- 
erage strength (s) and the length Ac- By changing the 
disorder strength we obtain behaviors such as those dis- 
played in Fig. 

Since interaction effects are negligible in the Gaussian 
limit, the characteristic length of density variations is 
set by the harmonic oscillator length a^, and not by 
the coherence length ^ as before, the latter being very 
large in that case. In this regime dominated by confine- 
ment, we observe that the shape of the density profile 
depends weakly on disorder in contrast to the Thomas- 
Fermi limit, for which this profile follows the variations 
of the disorder. This is particularly apparent in Figures 
|3Jc,rf), where disorder varies over a length scale smaller 
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FIG. 7: (a) Time evolution of the condensate density in the 
absence of disorder. p(0) is the stationary density and p{T) 
is the density after a time T = ^ . We have taken fj, = 2 and 
Said = 0.01 (b) Stationary profile of the condensate density 
in the presence of disorder. The corresponding values of the 
average disorder strength (s) , Ac and Ss are given in the inset. 
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FIG. 9; Left:(a,c,e) Plot of the disordered potential. The 
potentials in (a,e) vary on the same scale, whereas the po- 
tential in (c) varies on a smaller scale. The corresponding 
values of {s), Ac and 5s are respectively (0.14,80^,0.025), 
(0.25, 1.66a;,, 0.0846) and (0.31, 8a^, 0.105). Right: (b,d,f) 
Time evolution of the density corresponding to the potentials 
plotted on the left. The initial stationary value of the den- 
sity is denoted by puiiQ) and puiT) is computed at the time 
T = 25/uJz. We have used p — 2 and 2aid = 0.01. 



H<s> 



FIG. 8: Plot of the disorder averaged width (Luj) of the sta- 
tionary density profile in the Gaussian limit as a function of 
the average disorder strength p{s). The width is normalized 
by its value in the absence of disorder. For a given value of 
Ac, the width is averaged over 200 realizations of the poten- 
tial. The values of Ac are those used in Fig\^ and, in units of 
the harmonic oscillator length a,,, they range between 0.55az 
(lowest point) and Ba^ (highest point). We have used p — 2 
and 2a\d = 0.01. 

than the width of the density profile without leading to 
fluctuations of this profile. 

The density profile is well approximated by a off- 
centered Gaussian shape, 

p.(z) = Aexp(^-^i^i^^ , (32) 

for a large range of disordered potentials (see Figure 
13^6)). The amplitude A and the width are related to 
each other through the normalization. The average width 
{LJ) is a decreasing function of the disorder strength pi^s) 
defined by H17|l as represented in Fig|Sl Thus, the net ef- 



fect of disorder is to spatially localize the bosons inside 
a narrower Gaussian. 



B. Time dependent solutions 

The behavior of the stationary condensate density pro- 
file in the presence of disorder in the Gaussian limit dif- 
fers from the one obtained in the Thomas- Fermi limit. 
This difference shows up also in the time evolution of 
the density of the cloud after switching off the trapping 
potential. The short time expansion of the Thomas- 
Fermi cloud strongly depends on disorder, whereas in the 
Gaussian case, it does not. Moreover, in contrast to the 
Thomas-Fermi case, the zero point motion of the bosons 
is appreciable. The time evolution of the condensate den- 
sity after switching off the trap is presented in Fig. |^for 
different strengths of disorder. 

We first notice that on the same time scale, the density 
at the center of the cloud decreases more rapidly than for 
the Thomas-Fermi case (Figs|31and0J. This results from 
the non negligible kinetic energy of a Gaussian cloud and 
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[<s>, ^ , 8s] 
0.14,8a^,0.025 
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FIG. 10: (a) Time evolution of the spatial extension Lui{t) 
of the cloud for disorders of different strengths Vm and length 
scales of variation Ac- The disordered potentials are those used 
in Fig. \^and the corresponding parameters are indicated ({s), 
Ac, Ss) in the inset. Li^{t) is expressed in units of its value 
at t = 0. (h) Corresponding time evolution of the ratio TZ 
defined tn \31il . 



the weaker interaction between bosons. Figure 1101 dis- 
plays the time evolution of the average spatial extension 
{L^) of the cloud defined by (|28|l and the ratio of the 
average kinetic and interaction energies defined in ^\\\ . 
These two figures outline the difference between Thomas- 
Fermi and Gaussian time evolutions in the presence of 
disorder. The spatial extension in Fie tTUT a') does not 
show any saturation over comparable time scales, though 
it grows at a lesser rate with increasing the strength of 
disorder. Correlatively. the ratio TZ in Figure lTUT M grows 
at a much faster rate and it takes a longer time to sat- 
urate. We can summarize these observations by saying 
that though the cloud expansion is indeed prevented by 
the disorder potential in the Gaussian regime, the sup- 
pression is weaker than in the Thomas-Fermi regime and 
it happens on longer time scales. 



VI. SOLITON SOLUTIONS FOR AN 
ATTRACTIVE ONE-DIMENSIONAL 
BOSE-EINSTEIN CONDENSATE 

Having discussed the behavior of repulsive interacting 
bosons in the presence of disorder, we now turn to the 
case of an attractive solitonic condensate in similar sit- 
uations. As we shall see, the change of the nature of 
the interaction modifies the behavior of the soliton solu- 
tion with disorder as compared to the previous cases of 
Thomas-Fermi and Gaussian condensates. In contrast to 
Eq.® describing a repulsive interaction, Ea. (|ll|) involves 
one free parameter only {aid = As we have already 
mentioned, a change in aid only redefines the width of 
the soliton proportional to 1/^JJi. In what follows, the 
width is always kept less than ^. 



A. Stationary profiles 

We start with the study of the stationary solutions of 
Ea. (|ll|) with the addition of a random potential, namely. 



-^i<p{z) + dl<p{z) + Vd(l>{z) + 2|0|20 = 



(33) 



It is important to notice that, in contrast with previous 
cases, there is no trapping potential, so that in the ab- 
sence of disorder, the solution is invariant under trans- 
lations. Numerically, we start with a randomly chosen 
initial guess which, once iterated, gives a solution located 
around the initial trial function. The overall shape of the 
stationary solution turns out to be independent of disor- 
der, meaning that this shape can still be fitted with a 
hmction of the type AgSech{Bs{z — zq)), where and 



■ no disorder 
<s>=0.05, X^=^2i, 

. <s>=0.12, \=2.St, 

<s>=o.i9, >. =^2^, 




FIG. 11: Stationary density profile of a bright soliton in the 
presence of disorder. The chemical potential is fi = —20 and 
Old — —1. The disordered potential is characterized by \s\ 
(since jj. < 0) and Ac. < s > is the average of \s\. 

Bs — l/icj are respectively the amplitude and the in- 
verse width of the soliton. This feature appears clearly 
in Fig. Hll(l where the profile of the bright soliton has been 
plotted for several realizations of the potential. But, both 
the width and the amplitude depend on disorder as shown 
in Fig. H12|) which displays the behavior of the width for 
an increasing strength of disorder. We have also checked 
the dependence upon length scales Ac. Those features 
look similar to those obtained in the Gaussian limit. But 
they are essentially different. Whereas the soliton profile 
results from the comparison between kinetic and nega- 
tive interaction energies, the Gaussian profile is obtained 
from the comparison between kinetic and confinement 
energies. This difference will manifest itself in the time 
evolution of the solitonic condensate. 



B. Time-dependent solutions 

We now study the time evolution of the stationary solu- 
tions obtained previously, and not initial solutions given 
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FIG. 12: Disorder averaged width {Lu>) of the stationary pro- 
file of a soliton as a function of the strength of disorder. For 
one set of points (x), the strength Vm of the disordered po- 
tential is increased keeping fixed the spatial scale of variation 
Ac- For the other set of points (*), Ac is lowered which cor- 
responds to stronger fluctuations of the disorder, while Vm ts 
kept fixed. {L^) is expressed in units of its value in the ab- 
sence of disorder and, for each case, it is averaged over 200 
realizations of the potential. The average potential is charac- 
terized by \s\ since fj, < 0. The values of Xc are indicated in 
the figure. We have taken n = —20 and 2aid = — 1. 



by (|10|) unlike the case considered ini^. To this purpose, 
we first boost the soliton by giving it a finite (dimen- 
sionless) velocity Vs = 5. In the absence of disorder, the 
soliton travels a distance z = VgT over a time T without 
any change in its density profile. In the presence of a 
weak and smooth enough disorder, we observe that the 
soliton propagates retaining its initial (t = 0) shape, over 
distances comparable to the non disordered case. A weak 
disorder potential has thus a negligible effect on the soli- 
ton motion. For a stronger disorder strength (i.e., for 
a smaller value of Ac and a larger value of < s >), the 
time behavior is displayed in Figures IT^ a') and (&). In 
both cases, the soliton behaves classically and it becomes 
spatially localized, i.e. that it bounces back from high 
potential barriers typically higher than the kinetic en- 
ergy. However, we do not observe a significant change in 
the shape of the soliton. Its width fluctuates as the soli- 
ton travels through the disordered potential and bounces 
back and forth. When the strength of disorder is higher, 
the soliton motion is clearly not linear (Figure 6) ) . 
This kind of motion can be qualitatively explained by 
considering the soliton as a massive classical particle of 
mass mN , where m and N are respectively the mass 
and the number of atoms in the condensate. Deviations 
from the linear motion result from the spatially varying 
force exerted on the soliton by the disorder potential. 
This kind of description is valid as long as the disorder 
potential remains smooth over the width of the soliton. 



Similar behaviors have been discussed in the context of 
soliton chaos in spatially periodic potentials^, although 
the physical origin is different from the case discussed 
here. With the present stage of experiment^, such a be- 
havior can be verified by studying the time evolution of 
a bright soliton in an optical speckle pattern. 



la) 




FIG. 13: Time evolution of a boosted soliton in the presence 
of disorder. The chemical potential is fi — —20, the dimen- 
sionless velocity at t = is Va = 5 and aid ~ —1. The 
disorder potential is characterized by \s\ (since fi < 0) and 
Ac. (a) Fastly varying disorder with < s >= 0.12, Ss = 0.162 
and Ac — 2.8^. (b) Stronger but slowly varying disorder with 
< s >= 0.19, Ss = 0.372 and Ac = 12^. 



VII. CONCLUSION 

We have performed a detailed numerical investigation 
of stationary solutions and time evolution of one dimen- 
sional Bose-Einstein condensates in the presence of a ran- 
dom potential. Stationary solutions which correspond ei- 
ther to the attractive interaction bright soliton or to re- 
pulsive interaction Gaussian matter waves with repulsive 
interactions in the regime where confinement dominates, 
behave in a qualitatively similar way. In contrast, the 
stationary solutions that correspond to a repulsive inter- 
acting Thomas- Fermi condensate, depend strongly on the 
strength of disorder and on its spatial scale of variations. 

The time evolution of stationary solutions depends also 
significantly on the regime we consider. Although trans- 
port gets inhibited both for the attractive and repulsive 
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interaction, this occurs in a very different way. For the 
repulsive case the center and the edge of the cloud be- 
have differently and both arc ultimately localized in a 
deep enough potential well. In the interaction domi- 
nated Thomas-Fermi regime, the main part of the cloud 
remains localized and edges that correspond to low densi- 
ties and correlativcly weaker interactions, propagate fur- 
ther away. A study of the corresponding momentum dis- 
tribution of the cloud indicates a stronger localization of 
the matter wave in low momentum states for an increas- 
ing strength of the disorder potential. On the other hand, 
a moving bright soliton behaves very much like a single 
particle and it bounces back from a steep potential with 
its motion reversed. This behavior of a bright soliton 
may be contrasted against the behavior of a dark soliton 
in the presence of disorder which has been investigated 
recently^. 

For the values of the disorder strength and the non- 
linearity we have considered, we observe a behavior of 
solutions of the Gross-Pitacvskii equation that are mostly 
driven by the non-linearity, i.e., by interactions. Disor- 
der plays mostly the role of a landscape within which a 
classical solution evolves in time. We did not observe, 
for the relatively large range of disorder and interaction 
parameters we have considered, a behavior close to An- 
derson localization, namely where spatially localized so- 



lutions result from interference effects. Since disorder is 
expected to be stronger in one-dimensional systems, we 
may conclude that, for the currently accessible experi- 
mental situations, Anderson localization effects will not 
be observablai2iSSi^ due to the strength of the interac- 
tion term. Alternative setups are thus required in order 
to observe quantum localization of matter waves, having 
weak or zero interaction {e.g., by monitoring Feschbach 
resonances^i) . 

The signature of Anderson localization in the nonlin- 
ear transport of a BEG in a wave-guide geometry has 
been studied in2S. There, the transmission coefficient 
has been shown to be exponentially decreasing with the 
system size below a critical interaction strength. But the 
different types of disorder and of the matter wave den- 
sity at t = 0, make a direct comparison with these results 
difficult. 
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